Non-Markovianity as a probe for structural phase transitions in ion crystals 
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We investigate tiie back-flow of infoimation in a system with a second-order structural phase transition, 
namely a quasi one-dimensional Coulomb crystal. Using standard Ramsey interferometry which couples a 
target ion (the system) to the rest of the chain (a phononic environment), we study the non-Markovian character 
of the resulting open system dynamics. We study two different time-scales and show that the back-flow of 
information pinpoints both the phase transition and different dynamical features of the chain as it approaches 
criticality. We also establish an exact link between the back-flow of information and the Ramsey fringe visibility. 
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Introduction - It has recently been suggested that an open 
system approach [ 1 ] could be useful when studying properties 
of condensed matter physics, especially for systems present- 
ing critical features. In this approach the open system is re- 
garded as a quantum probe for the interacting many-body en- 
vironment. The idea, then, is to deduce some of the properties 
of the environment by monitoring the reduced dynamics of the 
probe only 10-01 • Here we apply this framework to Ramsey 
interferometry using a single spin in a Coulomb crystal [5- 
llt]. A Coulomb crystal can be depicted as a low-dimensional 
structure composed of several dozens of ions trapped in an 
anisotropic potentials with very large transverse confinement 
19]. This kind of atomic spatial aiTangement is routinely 
achievable in a standard Paul and Penning trap llTol - [l3ll . In- 
terestingly, Coulomb crystals exhibit slT -0 structural phase 
transition from a linear to a zig-zag equilibrium configuration. 
Such a transition can be tuned by a single experimentally con- 
trollable parameter, such as ion density or aspect ratio of the 
trap 1 14]. In this Article we present a method for accurately 
pinpointing the critical parameter determining the phase tran- 
sition. 

It has been shown that by time-monitoring the interfero- 
metric signal of the Ramsey protocol one can measure the 
autocoiTelation function of the crystal and detect the phase 
transition fl- Although this scheme is conceptually simple 
and experimentally feasible, two interesting questions are left 
open. If we look at the Ramsey protocol from an open sys- 
tem perspective, how does the dynamic of the single spin look 
like? Is there a way to signal the phase transition which would 
still rely on dynamical probing but, at the same time, be time- 
independent? Ramsey interferometry is made up of different 
unitary steps, making it impossible to define a unique Hamil- 
tonian describing the whole probe-chain interaction. Never- 
theless, one can imagine the single spin time-evolution as car- 
ried out inside a 'black box' and wonder about specific fea- 
tures of the consequent dynamical process. In this respect, 
we believe that a key feature addressing these questions is 
the (non-)Markovian character of the Ramsey interferometric 



scheme. 

Interaction between a single quantum system and its en- 
vironment leads to decoherence and dissipation, and the un- 
derstanding of such processes is of primary importance for 
both foundational and applicative reasons. The general the- 
ory of open quantum systems draws a line separating two 
major classes of processes: Maikovian and non-Markovian. 
In a Maikovian process the system irreversibly looses all its 
quantum properties to the environment. A rigorous charac- 
terization for this class exists, based on the ability of describ- 
ing the reduced system dynamics in terms of a semi-group 
of CPT dynamical maps |15]. A general definition of non- 
Markovian dynamics is more subtle and, as such, it has been 
subject of intense recent studies [16 - IS.] . Loosely speaking, 
non-Markovian behavior is characterized by memory effects, 
enabling the system to temporarily recover some of its gen- 
uine quantum properties. Here we focus on a novel definition 
for the non-Markovian character of an open system dynamics, 
originally introduced in Ref. [16]. The key idea to associate 
the distinguishability of two evolving system states to infor- 
mation flowing between the system and the environment: de- 
crease in the state distinguishability implies information flux 
from system to environment, while a reversed direction of 
information flow temporarily increases state distinguishabil- 
ity. This definition, and a subsequent measure quantifying 
the amount of non-Markovianity in a dynamical process, has 
been successfully applied to investigate memory effects in a 
number of different physical contexts, from photonic systems 
I.19i, .20,1 to cold atoms [3, .4J and spin-chains 

In this Article we investigate the dynamical map associated 
to Ramsey interferometry in a quasi- ID Coulomb crystal and 
study its non-Markovian behavior as the system approaches a 
phase transition. Our results show that the distinguishability- 
based measure of non-Markovianity is very sensitive to the 
critical point and, in this sense, can be used to probe the struc- 
tural phase transition of the Coulomb crystal and pinpoint 
the associated critical parameter unambiguously. Moreover, 
we find a direct link between this non-Markovianity measure 
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and the Ramsey fringe visibility, providing an experimentally 
feasible way for a novel characterization of this critical phe- 
nomenon. 

Coulomb crystal - A linear Coulomb chain is composed of 
ions with mass m and charge Q interacting via a repulsive 
Coulomb potential and confined inside a highly anisotropic 
trap. Let v and v, be the axial (x direction) and the transverse 
(y-zplane) trap frequencies. As shown in Ref. fT?], if we as- 
sume ring-type boundary conditions, a second-order structural 
phase transition form a linear to a planar zig-zag configuration 
takes place when 



vt < m \l (3) = vc. 



(1) 



where coq - -\/Q/ma^ and a is the equilibrium inter-ion 
distance. For v, < Vc the ionic equilibrium positions are 
n = (ia, 0, 0), with / = 1 , 2, . . . , A^. The ions oscillate har- 
monically around n resulting in small axial and transverse 
displacements from the equilibrium positions. In this spatial 
arrangement, the vibrational degrees of freedom of the chain 
are naturally uncoupled from each other. The normal modes 
of the chain can be labeled by the momentum k - 2jm/Na 
with n = 0, 1, . . .,N/2 and by their parity, v, is the largest 
transverse frequency of the normal modes excitation spectra. 
On the other side of the phase transition the ions reorganize in 
a zig-zag spatial configuration in the x-y plane where the new 
equilibrium positions are rj - (ia, (-iyb/2, 0). The transverse 
displacement b depends on both a and v, as shown in Ref. 
01411 . In this case the x and y vibrational degrees of freedom 
of the chain are coupled and the first Brillouin zone is now 
[0,7T/2a], with k = 27m/Na and /t = 0, 1, . ..,N/4. The exci- 
tation spectra, which displays multiple branches, are far more 
structured than the linear case. For a detailed derivation and 



description of the system see Ref. 111411 . 
Dynamics of the probe - In the following we concentrate 
on one specific ion of the chain, labelled j - I for conve- 
nience, and study its dynamics. Two internal states of this 
ion, {|e), \g)], form an open qubit system, and the rest of the 
ionic chain is taken to be its environment E. Since the relevant 
degrees of freedom of E are the vibrational ones, we can use 
the normal modes decomposition and picture £ as a bosonic 
environment at T = with a non-trivial discretized excitation 
spectrum. The qubit system S couples to the environment E 
via a resonant laser pulse that engineers the following interac- 
tion 
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where &^ - \e){g\,&^ - \g){e\, Q is the Rabi frequency of 
the laser, ul, ki are the laser frequency and wave-vector re- 
spectively, and yi is the ion position operator The general 
expression for yi can be found in Ref. [SJ. The target ion re- 
ceives a state-dependent mechanical 'kick' in the y direction, 
starts oscillating and excites all the transverse vibrational nor- 
mal modes of the chain. Operator e^''^'' ' is indeed a displace- 



ment operator for the transverse normal modes of the chain: 

^ik,h ^ (g) .) ^ (g) exp [afi] - a*bj) , (3) 
j j 

where ^j, bj are the creation and annihilation operators for the 
j-th normal mode of the chain, aj - irj ^iO(il2cL>jS \j, Lamb- 
Dicke parameter is rj - ki -\/h/m(L)o and matrix S ij, which 
realizes the normal modes decomposition, is defined in Ref. 
1 5]. The j index encodes all the quantum numbers defining 
the modes of the environment (momentum and parity). The 
Ramsey interferometric protocol we implement goes as fol- 
lows i^. Assuming oj mTl ^ 1, where is the largest 
eigenfrequency of the total system and Ti the laser pulse du- 
ration, we set Qr^ = 7r/4 in Eq.©, therefore applying a Jil2 
pulse, after which we let S and E evolve freely for a time t. 
The time-dependent dynamics is then governed by the follow- 
ing Hamiltonian 



(4) 



where ojj is the transverse frequency of the j-th mode. Fi- 
nally, we apply a -7t/2 pulse. The total evolution operator 
reads as follows 



U(t) = UiMT(-n/2)Uo{t)UiNT{7r/2). 



(5) 



It is important to remark that t is the time elapsed between the 
two pluses. Using the Bloch-sphere representation [6, (p) for 
S , we choose a generic pure initial state of the form 



l<A/> 



cos(^)|e)+e'Mn(^)|g)](g)|0>, (6) 



where |0) is the total phononic vacuum state, which is defined 
as |0) - \0j}. After applying ^ to we obtain the 
following final system-environment state 

UitMi) = 1^/) = i [\g,Xg(t)} + \e,Xe(t)}] , (7) 



where 



\Xg} = cg\0} + ic,\ -a) + CgD(-a)\a(t)} - ic,\ - a(t)} 
\Xe} = q|0) + icg\a) + CeD(a)\ - a{t)) - icg\a(t)) 



(8) 



and we have introduced the short-hand notation c^ - 
cos {912) , Cg = e'"^ sin (61/2), \a) = (g)^. \aj), aj(t) = aje-''^^' 
and D{a) = (^jD(aj). Since we are interested in the re- 
duced system dynamics only we trace out the environment 
degrees of freedom to get the state of the probe qubit ps (r) = 
Peeit)\e}{e\ + pgg(t)\g){g\ +p,g(t)\e){g\+ h.c. 

The overall process admits two complementary viewpoints: 
If we set the initial time to = -Ti, then the system time- 
evolution is dictated by a complex quantum map, to which 
we cannot assign a unique global Hamiltonian describing its 
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action. If, instead, we set fo = we focus on the free evo- 
lution of an initially correlated system-environment state. Ei- 
ther way, we are effectively engineering and simulating a non- 
trivial open system dynamics where both dissipation and de- 
coherence processes can take place. In the following we will 
quantify the degree of non-Markovinaity associated to this dy- 
namical process. 

Non-Markovianity and criticality - We start by recalling the 
non-Markovianity measure AT introduced in Ref. [16]. Given 
two quantum states p\,p2 the distinguishability between them 
is quantified by the trace-distance D - tr|pi -p2l/2. A 
generic Markovian dynamical map pit) - Ktp{to) is contrac- 
tive, that is, it continuously decreases the trance distance be- 
tween two initially distinguishable quantum states. If at least 
two initial states exist such that this condition is violated for 
some time intervals then the map is non-Markovian. The non- 
Markovianity of a dynamical process described by a map A, 
can be quantified by the following integral 



7V(A) 



: max 

Pl,2(0) 



dto-(t,pi^2(0)). 



(9) 



where cr(f,pi,2(0)) - j,D(p\{f),p2(t)) is the rate of change of 
the trace distance, that is, the information flux.. As shown in 
Ref. 1231], the maximization in Eq.(|9]l can be restricted to pairs 
of initially pure orthogonal states pi(0), p2(0) for the case of 
a two-level system. These states, which lie on the boundary 
of the state space, are represented by antipodal points of the 
Bloch-sphere surface. This also explains the ansatz (|6]l for the 
initial state of the probe. 

We now study the behavior of AT as a function of the tuning 
parameter A - v,/vc - 1. The dynamical map for the reduced 
density matrix reads as follows 

ps(t) = A,ps(0) = tTE [UiDps^EUHt)] , (10) 

where U{t) is given in (|5]l. Given its complexity, it is impos- 
sible to perform the maximization in Eq.(|9]l exactly when the 
number of ions is large. However, we have shown numeri- 
cally that the maximizing pair is formed by the eigenstates of 
cr,, which we label |-). The corresponding trace distance 
reads as Dopt(t) - \p%\t)-p'eg\t)\- Interestingly, the dynamics 
of |+), |-) is purely dephasing, i.e., the optimal pair does not 
exchange excitations with the environment. 

In the following we show that the dynamics of the trace 
distance displays a sensitivity to the tuning parameter A: the 
closer to criticality the chain is, the more damped the time- 
oscillations of D are. Its temporal behavior is very similar to 
that of the visibility of the Ramsey interferometric signal jst] 



n^(f) = exp[-A(f)], 



(11) 



where A(f) = 2 Yjj Iq'jP sin^ (a>jt/2y Indeed, the trace distance 
can be formulated in terms of the visibility 'V(f): 



Dopt(t) 



1 + 2 cos [Bit)] l^^it) - ^ j + V(f) + 2^4 ^ 

(12) 
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FIG. 1. Short-time-scale dynamics, (a) Time-evolution of the trace 
distance for the maximizing pair {|+>,|->) for N = 200, far from 
the critical point (black line, A = 0.1) and very close to it (red line, 
A = lO"'*). (b) Non-Markovianity measure N for short time-scale 
truncation as a function of A for A' = 100 (blue solid line) and N = 
1000 (green dashed line). 



where ^ = g-S/lfjI'/z ^jj^j g^^f-^ _ Yjj\aj\^ sinicjjt). This con- 
nection is one of our main results, establishing an exact ana- 
lytical link between the experimentally measurable fringe vis- 
ibility 'Vit) and the non-Markovianity N of the probe qubit. It 
is worth noting that the term 'Vit) - ^'^ j'Vit) accounts for the 
overlap between time-dependent coherent states of opposite 
sign and the initially laser-generated coherent states. This fea- 
ture is a consequence of the optimal pair of states and shows 
notable similarity to the critical model of Ref. [|2j]. Let us 
remark that the above formula is general and specific of the 
dynamical steps (|5]) only. The nature of the chain, that is ei- 
ther linear or zig-zag, is fully encoded in ^, A(f) and Bit) only. 

Short time-scale - We first analyze the short time behav- 
ior of the trace distance. In an experiment-based spirit, we 
set an upper bound to the elapsed time, which roughly 
corresponds to 250/is (this time will also be the upper inte- 
gration limit in (|9j) [10], define t = wyf and look at the 
time-evolution of the trace distance of the maximizing pair, 
FigdJa). The behavior of A/' as a function of A is shown in 
FigEb) for = 10^, 10^. Coordinate A = represents the 
critical point and the maximum distance from it that we con- 
sider is 10"', whereas the minimum is 10"'^. There is a clear 
and rather abrupt change in the behavior of N at the critical 
point, which coincides with the appearance of the structural 
phase transition. Values of N are different on different sides 
of the critical point. However, approaching criticality both 
Na^q- and Na^o* converge to the same value, a non-zero ab- 
solute minimum. This result fully reflects the time-evolution 
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FIG. 2. Long-time-scale dynamics, (a): time-evolution of the trace 
distance for the maximizing pair [\+), |-)), for N = 300, far from the 
critical point (black, A = 0.1) and very close to it (red, A = 10"^). 
(b): Non-Markovianity measure N for long time-scale truncation as 
a function of A for A' = 300. 



of the trace distance, as shown in FiglHa). As the environ- 
ment approaches the critical point the dynamics of the probe 
becomes more and more damped, resulting in a Markovian be- 
havior. Interestingly, the behavior shown in FigJTJb) is char- 
acteristic of a second order phase transition. Indeed, even if N 
is a continuous function when A — > O'', its derivative is clearly 
not. A similar feature is found when one studies derivatives 
of thermodynamical and statistical quantities in presence of a 
classical and quantum phase transition. 

Long time-scale - The next step is the investigation of the 
long-time scale behavior, corresponding to fM ~ 5 ms. The 
dynamics of the trace distance in this case is displayed in Fig. 
|2{a). Contrary to what happens in the short-time-scale regime, 
the closer A to zero, the larger the amplitude of the slowly 
damped revival peaks in D. Consequently, we expect the non- 
Markovianity measure to increase as the chain approach the 
critical point. Fig|2] fully confirms this prediction as N dis- 
plays a cusp-like maximum when A — > O''. At criticality the 
mechanical instability of the chain is driven by the so-called 
soft mode [14]. This mode is associated to the smallest fre- 
quency possible in the environment spectrum, which corre- 
sponds to k = 71 1 a and vanishes as A — > 0=^. Hence, close to 
the phase transition and for long enough times, the coupling 
between this almost stationary mode and the single spin will 
dominate over all the others, leading to a weakly damped os- 
cillatory dynamics of the probe. The action of the phononic 
background, which would force the system to lose informa- 
tion at short-time-scales as a consequence of the instability, 
is here overruled by an effective one-to-one coupling between 
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FIG. 3. Short time-scale non-Markovianity measure Ncr at criticality 
as a function of the number of ions A'. 



the target ion and the soft mode. It is important to remark that 
the short and long time-scale regimes are totally compatible, 
with the common feature being the strong coupling between 
the probe and the soft mode at criticality. The main differ- 
ence lies in the fact that for short times the probe is able to 
resolve only the high-energy part of the environment spec- 
trum. On the contrary, if we wait sufficiently long time the 
revivals in the dynamics of the trace distance, due to the cou- 
pling with the soft mode, will appear leading to a the peak in 
non-Markovianity. 

Finite size effects - Finally we address how finite size ef- 
fects affect our results. To this aim we study how the short- 
time scale non-fVIarkovianity measure at criticality scales with 
the number of ions in the chain. The results are reported in 
Fig. |3] Ncr saturates to a non-zero value very soon. Finite-size 
effects, which lead to a larger value of the non-Markovianity 
measure, are relevant only for relatively small and no ap- 
preciable variation of Na- is detectable for > 100. This is 
an indication that even in thermodynamical limit, the flow of 
information from the probe qubit to the rest of the chain is 
not complete. Loosely speaking, a minimum amount of the 
information initially stored in the probe remains trapped and 
does not get lost ever. As previously pointed out, for short 
times only the high-energy part of the spectrum (small k) con- 
tributes to the dynamics of the probe. This portion of the en- 
vironment spectrum is essentially flat {d(jL>{k) j dk\k=Q - 0) and 
so all these modes dephase and rephase almost in sync; the 
system only leaks information. However, this happens on av- 
erage only as the dephasing and rephasing cycles repeat many 
times within the short time interval we consider. This, in turn, 
implies that some of this information cannot flow out. It could 
do so completely if there were not the short-time-scale recur- 
rences. These are associated to mechanical excitations going 
back to the system, resulting in revival peaks in the dynamics 
of the trace distance. The typical time for these revivals was 
estimated as t » 140. 

Conclusions - We have shown that the back-flow of 
information in single-spin Ramsey interferometry is a witness 
for the linear-to-zig-zag structural phase transition of a quasi 



5 



1-D Coulomb crystal. Abrupt changes in the behavior of 
TV are observed when v, — > +Vc, leading to extrema which 
clearly pin-point critic ality. The nature of these extrema 
consistently reflects the dynamical features of the chain, 
which are responsible for driving the phase transition. Again, 
we want to stress that the above analysis is experimentally 
feasible. By time-monitoring the visibility of the Ramsey 
fringes, one can compute the non-Markovianity N of the 
interferometric scheme using Eq. rtTZt . This method relies on 
dynamical probing and, at the same time, characterizes the 
phase transition in a time-independent way. 
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